Electroneutrality and the Friedel sum rule in a Luttinger liquid 
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Screening in one-dimensional metals is studied for arbitrary electron-electron interactions. It is shown that for finite-range 
interactions (Luttinger liquid) electroneutrality is violated. This apparent inconsistency can be traced to the presence of external 
screening gates responsible for the effectively short-ranged Coulomb interactions. We also draw attention to the breakdown of 
linear screening for wavevectors close to 2kF. 
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Screening is one of the most important and useful 
concepts in condensed-matter physics [0,^. If some ex- 
ternal charge is brought into a conductor, the internal 
charge carriers will reorganize with the new distribu- 
tion of charge eliminating the electric field at large dis- 
tances. The screened potential set up by the external 
charge together with its screening cloud is then rather 
short-ranged. Typically, one ends up with only weakly 
correlated systems, despite of originally long-range elec- 
trostatic forces. 

Most theories employ linear screening as working hy- 
pothesis, where the effects of (possibly time-dependent) 
external test charges onto the conduction electrons are 
determined by linear response theory. Given the valid- 
ity of linear screening, the wavevector- and frequency- 
dependent dielectric function e{q,uj) contains all relevant 
information about screening. An important result of the 
theory is the Friedel sum rule which states that the 
total electronic screening charge exactly compensates any 
external (impurity) charge brought into the system. This 
charge neutrality requirement on large scales arises be- 
cause in equilibrium there can be no net electric field at 
large distances. The validity of the Friedel sum rule is 
usually taken for granted, generally by referring to the 
analysis in Ref. Q where this was proven explicitly for 
the Anderson model. 

In this paper, we discuss screening and the Friedel sum 
rule for interacting electrons in one dimension (ID). At 
low energy scales, provided no lattice or spin instabil- 
ities are present, the properties of ID fermions can be 
described by the Luttinger liquid model Q (or slight gen- 
eralizations thereof, see below). The Luttinger liquid is 
a strongly correlated ID metal which does not support 
the existence of Landau quasi-particles. It is of impor- 
tance for a number of applications of current interest, 
e.g., quantum wires in semiconductor heterostructures in 
the limit of one transport channel Q, transport in car- 
bon nanotubes [Q, or quasi- ID organic conductors [||, 
to mention a few. Here we show that in a Luttinger liq- 
uid the screening charge does not balance the impurity 
charge. Nevertheless, Friedel's phase shift sum rule 0] 



remains valid in terms of an adequately defined phase 
shift. 

Our investigation of screening in ID is based on the 
standard bosonization method |^]. From a comparison 
with alternative techniques, this method is known to give 
a proper description of ID fermions at low energy scales 
(we consider zero temperature below). Since spin and 
charge are decoupled in a Luttinger liquid, it is sufficient 
to study only the spinless case in the following, with the 
same conclusions applying to spin-^ electrons. The low- 
lying excitations in a spinless Luttinger liquid arc de- 
scribed by a bosonic phase field 0{x), in terms of which 
the electron density operator can be written in the form 

p{x) = — + -^dxO{x) + — cos[2kFX + Vl7r0{x)] . 

(1) 

The first term describes the mean charge density (which 
is supposedly neutralized by a positive background) , the 
second term gives long-wavelength [q ~ 0) fluctuations, 
and the last term yields rapidly oscillating {\q\ ~ 2/cf) 
contributions. Putting Ti — \, the Luttinger liquid Hamil- 
tonian is then given by ||,^ 
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(2) 



+ — / dxdx' dJ{x)U{x - x')d^,e{x') , 

ZTT 



where n(x) is the canonically conjugate momentum to 
9{x) and vp the Fermi velocity. The potential U{x) 
can describe either an unscreened Coulomb interaction, 
U{x) ~ or an externally screened finite-range po- 

tential which arises due to the presence of mobile charge 
carriers close to the ID metal, e.g., on screening gates or 
other nearby chains. To simplify notation, we shall make 
the inessential assumption that U{x) is sufficiently long- 
ranged such that its Fourier transform U{ci) has a very 
small component at q = 2k p. Then electron-electron 
backscattering can be neglected, as implied in Eq. (||). 
Even if this should not be the case, the bosonization 
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technique can still be applied and yields qualitatively the 
same results. The Luttinger liquid model in the strict 
sense is obtained by effectively using a local interaction 
U{x) = U{0)6{x). Here we shall employ the usual dimen- 
sionless Luttinger liquid interaction parameter g defined 
as 



5 = [1 + U{0)/ttvf] 



^1/2 



(3) 



for all finite-range interactions. Then g = 1 is the non- 
interacting limit, while for repulsive interactions we have 
g < 1. In the absence of screening gates, g approaches 
zero in an infinitely long system [^J9| . 

To study screening properties, we now consider some 
external time-dependent charge distribution eQ{x,t) 
brought into the system. The interaction with the ID 
metal reads 



Hqit) = / dxdx' Q{x,t)U{x~ x')p{x') 



(4) 



In view of the representation (|l|) for the electronic den- 
sity, there are two contributions. The first comes from 
the g ~ component, and the second from the q ~ 2kp 
part. We note that the interaction potentials in Eqs. (|^) 
and are the same because we deal with the internally 
unscreened, microscopic interaction at this stage. 

Let us first discuss the long- wavelength (|g| ^ 2k p) re- 
sponse of the electrons. Ignoring the 2kF part in Eq. (|^), 
the now Gaussian Hamiltonian yields straightforwardly 



{piQ,^)) = 2 — 2r^Q[i,i^) 

-K bJ'^ — OJ (q) 
with the plasmon dispersion relation 



(5) 



uj{q) = VF\q\yi + U{q)/TTVF ■ 

Apparently, in the long- wavelength limit, the Luttinger 
liquid model implies linear screening 

{pii,^)) = U{q)x{q,i^)Q{q,i^) 

with the polarizability x(<Z, The response of the elec- 
trons to Q{x, t) is thus fully described by a dielectric 
function. One finds from Eq. (||) and the definition |Q] 



e-\q,u:) = l + C/(g)x(g,w) 
the small-g result 

e (g,w) = l + — 

TT cj^ — U)'^{q) 

In the static case, = 0, this yields 
e{q) = 1 + U{q)/TTVF ■ 



(6) 



(7) 



We mention in passing that for large impurity charge 
eQ the bosonization approach breaks down ]lO| . For in- 
stance, Eq. (||) would incorrectly predict that the electron 
density becomes negative for sufficiently large Q. 

One can now define the internally screened interaction 
potential Ucs{q) = U{q)/e{q) which determines the effec- 
tive potential between two charges J^. From Eq. (^, its 
long-wavelength part is Ucs{q) = U{q)/[1 -\- U{q)/TTVF], 
which gives for a finite-range interaction C/eff (g) — g^U{q) 
as g 0. For an externally unscreened l/\x\ interaction, 
one has U{q) — 2e^| ln(g(i)|, where d is the width of the 
ID channel Apart from a hard core at small dis- 
tances, this leads to the large-distance behavior valid at 
|a;| > d 



|x|ln|a:/d| 



(8) 



Therefore the long-range character of the interaction is 
not significantly reduced. The only logarithmic suppres- 
sion of the l/\x\ law explicitly demonstrates the very 
weak screening in ID. 

The condition for perfect screening 



-i(g->0,c^-0) 







is seen to be violated in any finite-range model. This 
follows directly from Eq. (j^) since e(g ^ 0) = The 
implications are best discussed for a point charge sitting 
at a; = 0, i.e., Q{x,t) — QS{x). The corresponding long- 
wavelength response is given in Eq. (^). For the total 
screening charge, eQs — e J dx{p{x)) , this leads to the 
strikingly simple result 



Qs 



(9) 



where g has been defined in Eq. (||). We stress that this 
relation holds for any finite-range Coulomb interaction. 
Asserting that the 2kF Friedel oscillation in the charge 
density does not contribute to the total screening charge 
(see below), we observe that only a fraction 1 — g^ < 1 
of the external charge Q is screened by the conduction 
electrons. Therefore the electroneutrality condition for 
impurity plus screening charge, Qs + Q = 0, is appar- 
ently violated in models with a finite-range interaction. 
Of course, this reasoning carries over to lattice models 
with effectively short-ranged interactions, e.g., the ID 
Hubbard model. For a long-range l/\x\ interaction, the 
parameter g effectively approaches zero, and electroneu- 
trality is then seen to hold. 

The result can also be obtained by a phase shift 
consideration. Forward scattering due to a point-like im- 
purity charge Q5(x) in Eq. (|4|) can be eliminated by the 
standard unitary transformation 

U = exp{— i\/7r / (ia;a(x)0(x)} , 
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where (j){x) is the dual field to 9{x) 
transform of a{x) is 



a{q) = - 



U{q)/TTVF 



1 



Q 



and the Fourier 



(10) 



U{q)lTiVF 

Comparing with Eq. (|^) for lj — the induced electronic 
density is simply {p{x)) = a{x). The unitary transfor- 
mation U now leads to a phase shift appearing, e.g., in 
the 2kF part of Eq. (|l|), which takes the form 

ri{x) = vr / dx' sgn(x — x')a{x') . 



Defining the asymptotic phase shift rip — rj[x ^ go), we 
find 



= / dxa{x) = tif/t^ 



(11) 



Despite the apparent violation of electroneutrality, 
Friedel's phase shift sum rule |^ is seen to hold. Clearly, 
the phase shift -qp characterizes some screened impurity 
charge and not the bare charge Q brought into the sys- 
tem. Finally, using Eqs. (^) and (^), one may verify 
Eq. (|) again. 

We mention that the conventional Fermi liquid case is 
not directly included in Eq. (^ as the simple limit g — > 1. 
For a Fermi liquid, one assumes that quasi-particles with 
good screening exist and then adds a local potential scat- 
terer in order to derive the Friedel sum rule Its scat- 
tering strength is related to a phase shift rjF, and the 
screening charge is Qg = t{fIt^ as in Eq. (|ll|). By inter- 
preting rjF /tt as the impurity charge, the Friedel sum rule 
is then in fact imposed as a consistency relation ensur- 
ing electroneutrality of the system. In contrast, putting 
(7 = 1 for the Luttinger liquid model would imply a non- 
interacting system (Fermi gas) rather than a Fermi liquid. 

The physical reason for the apparent failure of elec- 
troneutrality in finite-range models are induced charges 
outside the ID system, e.g., on external screening gates, 
which cause the finite range of the interaction. These 
other conductors also contribute to the total screening 
charge. To give a concrete example, consider the gated 
ID quantum wire shown in Fig. 1. For a wire of width d, 
the presence of a two-dimensional gate at a distance D 
leads to a short-range interaction characterized by 



2p2 

g= ll + ^ln(2L>/d) 
ttvf 



-1/2 



The induced ID charge density [integrated over the y- 
direction] on the gate, epcix), is obtained as 



D f 

Pg{x) = -— / dx' 



q{x') 



+ {x- x'y 



where q{x') is the total density in the wire (including the 
impurity) . Integration over x gives straightforwardly the 
total induced charge on the gate 



Qg = j dx pg{x) = -{Q + Qs) ■ 
In effect, the electroneutrality condition in the form 

Qs + Q + Qg = Q (12) 

is then restored for the total system including the screen- 
ing gates. Ignoring the screening gates implicitly used to 
derive the Luttinger liquid model is thus responsible for 
the modified condition (||) within the ID system. Paren- 
thetically, we note that if the charge Q is not put directly 
into the ID system but some distance away, the screening 
charge Qs is not given by Eq. (^ anymore, yet Eq. ( p^ ) 
will still hold. 




FIG. 1. Schematic view of a ID quantum wire with 
sliort-range Goulomb interaction due to the presence of a 2D 
screening gate located a distance D away from the wire. The 
bare impurity charge is Q, the direct screening charge within 
the ID system is Qs = —(1 — g^)Q [see Eq. (|^], and the 
induced charge on the screening gate is Qg- 

The underscreening of an external charge brought into 
the ID system should also be experimentally observable. 
If charge is injected into, e.g., a carbon nanotube con- 
stituting a perfect experimental realization of a ID con- 
ductor [01, the resulting image charge on nearby exter- 
nal screening gates can be detected by capacitance spec- 
troscopy [jll| or by using highly sensitive single-electron 
transistor (SET) electrometers on top of a scanning probe 

tip g. 

So far we have discussed the long-wavelength part of 
the electronic response only. Turning now to the 2kF 
part [ p^ and assuming linear screening, we have to com- 
pute the polarizability which is essentially the double- 
Fourier transformed density-density correlation function 
of the unperturbed conductor 



X((7, it^) = - 



*«-(T,p(.T,T)p(0,0)) . (13) 



Here Tr is the time-ordering operator in Euclidean time, 
and one has to analytically continue Eq. ( p^ ) to real fre- 
quencies, iio — > cj + iO+, in order to obtain x('?, ^) needed 
in Eq. (^. The 2kF part of x for a Luttinger liquid is 
found to read 

2g-2 

x{q,i^) = -— V 

■kvf 

p.s=± 



ISLO 

vpkp 



kf 



^2p 
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isLU — vplq — 2pkF\ 
isuj + VF\q — 2pkp\ 

where F denotes the hypergeometric function and Cg — 
4-sr(2 - 2g)/r{g)T{2 ~ g). In the static case, this gives 
algebraic singularities 



X{q) = E 



hp 



2g-2 



(14) 



with the numerical constant 



Co 



V^r(2 - 2g) 
2T{g)V{i/2-g) 



From these algebraic singularities one would infer a 
Friedel oscillation decaying as (p(x)) ^ cos{2kpx)x^^'^^ 
and a similar contribution to the screened interaction po- 
tential Ues{x). However, this represents only the first 
order in the perturbation expansion for the Friedel os- 
cillation and determines merely the short-distance be- 
havior, while the long-distance behavior of the Friedel 
oscillation necessitates a calculation in all orders of the 
impurity strength ]T^ , p^ . An important implication is 
the breakdown of linear screening for the 2kp electronic 
response. This breakdown occurs for arbitrarily small 
impurity charge eQ at wavevectors close enough to 2k p. 
Following the results of Ref. ||lj], the singularity expo- 
nent 2g — 2 ioY q ^ 2kp in Eq. ( p^ ) is turned into 
g — 1. As a consequence, the effective screened poten- 
tial as well as the Friedel oscillation asymptotically decay 
as ^ cos(2A;j;-a;)a;~^. Due to the intrinsically nonlinear 
screening, the dielectric function is of rather limited use 
for wavevectors close to 2k p. 

The bosonization approach naturally separates the 
density operator (p, and therefore also the electronic 
screening response, into a slow and a fast 2kp part. The 
total screening charge is determined by the q — Q compo- 
nent of the induced charge density, which in turn is exclu- 
sively given by the slow part (|^) . Therefore the Friedel os- 
cillation obtained from the bosonized 2kp part of Eq. (|l|) 
does not contribute to the total screening charge. In a 
microscopic calculation, one will in general not be able to 
separate the slow and the fast components so nicely, but 
within the bosonization approach, a quite simple deriva- 
tion of the total screening charge (||) is possible. 

To conclude, we have investigated screening in one di- 
mension. We have shown that electroneutrality is not 
obeyed in models with a finite-range (screened) Coulomb 
interaction. In a ID metal, the total induced screening 
charge is given hy Qs ~ —{l~g^)Q, where Q is the 
impurity charge and g the Luttinger liquid interaction 
parameter. To resolve this apparent inconsistency, one 
needs to take into account induced charges on screening 
gates. 
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